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We present a mean-field calculation of the phase diagram of a simple model of 
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Q\ . localized moments, in the hexagonal uranium heavy-fermion compounds. The model 

o , 

. considers a non-Kramers quadrupolar doublet ground state magnetically coupled 



with a singlet excited-state, favoring in-plane van-Vleck magnetism, as has been 
conjectured for UPta. The Hamiltonian which defines the model is Heisenberg like in 



C , both, magnetic and quadrupolar moments. No Kondo effect physics is included in the 



calculations. Among our main results are: (i) for zero intersite quadrupolar coupling, 



^ , the magnetic order is achieved by a first order transition above a critical intersite 

magnetic coupling value which becomes second order at higher coupling strengths, 
(ii) for finite intersite quadrupolar coupling, at temperatures below a second order 
quadrupolar ordering transition, the minimal magnetic coupling value is increased, 
but (a) the magnetic ordering temperature is enhanced above this value, and (b) the 
ordering of first and second order transitions in the phase diagram is reversed. By 
considering the general structure of the Ginsburg-Landau free energy, we argue that 
the Kondo effect will not modify the shape of the phase diagram, but will modify 
the quantitative values at which transitions occur. 
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I. INTRODUCTION 



The unusual commensurate, weak moment magnetic order present in heavy electron 
superconductors such as UPts [and its related alloys U(Pti_a;,Pda;)3 and (Ui_2:)Th2;)Pt3] and 
URu2Si2 has so far resisted explanation in any simple terms. It is difficult to believe the 
small moment order represents itinerant spin density wave formation, since, for example, 
UPts displays no appropriate nesting features on the Fermi surface at the appropriate anti- 
ferromagnetic wave vector. It is difficult to account for the order in terms of frustration 
effects since the ordered structures are remarkably simple. In the case of URu2Si2, moreover, 
there is evidence fl]] that the magnetic order parameter is not the 'primary' one, in the sense 
that despite the large entropy involved in the 17. 5K magnetic transition (the entropy involved 
is larger than in the 1.2K superconducting transition) the effective moment is too small to 
be consistent with the size of the ordered moment. Indeed, it has been suggested that 
some other "hidden" order drives the transition (such as quadrupoles) . 

Recently results coming from different directions altogether make it prudent to consider 
quadrupolar (and octupolar) ordering possibilities in these materials: it has been noted that 
the quadrupolar Kondo model beheved to be previously restricted to cubic uranium based 
materials is in fact possible for hexagonal (UPts) and tetragonal (URu2Si2) materials as well 
0]. In this model the heavy fermions arise from local uranium quadrupole moments associ- 
ated with ground crystal field doublets quenched by (predominantly) the orbital motion of 
itinerant electrons. Thus, the low lying degrees of freedom may be essentially quadrupolar. 

The relevance of the quadrupolar Kondo model to this class of materials has been put on 
ffimer footing by the discovery of the cubic material Yi.^^U^Pds, which appears to display 
a dilute limit quadrupolar Kondo effect [Q. Moreover, all known heavy fermion supercon- 
ductors possess the appropriate symmetry groups to yield a quadrupolar Kondo effect or 
magnetic two-channel Kondo effect (CeCu2Si2), a point bolstered by the recent discovery of 
two hexagonal heavy-electron superconductors (UPd2Al3 and UNi2Al3) which also possess 
magnetic order 0. 
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In the quadrupolar Kondo picture, a context for weak moment magnetic order can arise 
naturally from two distinct sources: (i) for cubic symmetry, all magnetic character is in 
excited states, and for hexagonal and tetragonal symmetry, in plane magnetic character is 
in excited states. Thus, the magnetic response is van Vleck in character yielding a built in 
stability against magnetic ordering and a region of weak moment order, (ii) For hexago- 
nal and tetragonal symmetry, the c-axis character of the local uranium doublets is weakly 
magnetic, but primarily octupolar. 

In this paper, we shall focus our attention on (i), with the particular goal of examining 
how quadrupolar order can induce or at least enhance magnetic ordering tendencies. While it 
is well known that quadrupolar order will always be induced as a secondary order parameter 
at a magnetic transition, it seems to be less well known, at least in the heavy fermion field, for 
magnetic order to be induced or enhanced by the coupling to quadrupolar degrees of freedom. 
Such physics is known in UO2, for example, where the antiferromagnetic phase transition is 
first order due to a coupling between magnetic and quadrupolar order parameters 0. 

We have performed mean field calculations on a simple model of potential relevance to 
the hexagonal uranium heavy fermion superconductors. In this model, each uranium ion 
possesses a doublet quadrupolar ground state magnetically coupled to an excited singlet 
state. The uranium ions are placed on a lattice with the structure of UPts (though the 
point symmetry is taken to be D^h rather than -D3/1). We then turn on nearest neighbor 
in-plane and out-of-plane quadrupolar and magnetic coupling and study the resulting phase 
diagram. We find a rich structure, ranging from van Vleck magnetism with a tricritical 
point and critical exchange strength to alternately suppressed and enhanced magnetic order 
below a second order quadrupolar transition. We discuss the phase diagram using both the 
full mean field calculations and Landau theory in regions of validity. 

While the model does not produce a full explanation of the ordering in bulk UPta (for 
which the specific heat shows no transition and the neutron elastic scattering peaks are not 
true resolution limited Bragg peaks) it may be of relevance to the Pd and Th doped alloys 
and to the new hexagonal superconductors. In any case, it represents the beginning of an al- 
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ternative view of collective phenomena in the heavy fermion materials which may ultimately 
provide a unifying view of the heavy fermion state, superconductivity, and magnetism. 

The outline of the paper is as follows: Section II provides an overview of the model. 
Section III develops the mean field approximation to this model. Section IV discusses our 
results in terms of Landau theory, paying special attention to the role of the orientation 
dependent coupling between quadrupolar and magnetic order parameters. Finally, in Section 
V we conclude and point out explicit experimentally relevant features of our model. 

II. MODEL 

In the section II-A we write a single-site Hamiltonian, coupling the quadrupolar moment 
of the ground state of the uranium-ion with the electric-field gradient and the magnetic 
moment with an applied magnetic field. In the next section we introduce intersite coupling 
of the XY-model form for both quadrupolar and magnetic degrees of freedom. Within the 
mean-field approximation that Hamiltonian is equivalent to that of section II-A, but with 
coupling coefficients which depend self-consistently on the mean magnetic and quadrupolar 
moments. 

A. Single-Site Hamiltonian 

Having in mind the motivations given in the Introduction, we define a single-site model 
in which the uranium ion has a non-Kramers quadrupolar doublet ground state magnetically 
coupled only with an excited state. Such a crystal-field scheme may arise for a U^"'", 5/^, 
J = 4 at a site of hexagonal symmetry. Fig. defines the level structure of the model, where 
the doublet \E±) and the singlet \B) are defined by 

\E±)=u\±A)+v\t2) , (1) 

and 
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m ^ . (2) 

where \m) are the eigenstates of the z-component of the total angular momentum J. u and 
V are constants. The relevant matrix elements for the model we are defining are 

{B\J±\E^) = 2u + V7v = a, (3) 

and 

{E±\Jl\E^) = 2V7v = b. (4) 

Since + = 1, the constants a and b defined here are not independent. 

The other J = 4 multiplet levels are not considered, so the model will be restricted 
to temperatures less than or equal to the first excitation energy A. It is clear from the 
level structure considered that this model favor in-plane van Vleck magnetism, and c-axis 
Pauli magnetism, as has been conjectured for UPts 0. The motivation for focussing on the 
doublet level is that it is the only degenerate level in hexagonal symmetry for ions, and 
we desire the internal degrees of freedom associated with the doublet level to allow for the 
Kondo effect. The formation of the many body Abriksov-Suhl-Kondo resonance is the basic 
paradigm for the source of heavy fermions in these materials. 

The doublet ground state corresponds to the E2 representation of the point-group D^h, 
where the eigenf unctions transform like — y"^ and xy. This doublet has magnetic char- 
acter along the c-axis, but only in-plane quadrupolar tensors may couple the doublet levels 
together. This implies we can have a quadrupolar Kondo effect. Using the Stevens op- 
erator representation [^, the corresponding quadrupolar operators are (proportional to) 
Qxx = ~Qyy = Jl ^ ^-iid Qxy = Qyx = {JxJy + JyJx) l'^- The magnetic moment operator 
is M = gjfisJ, where gj = 4/5 and fis is the Bohr magneton. The doublet is connected by 
Jx, and Jy to the excited \B) singlet. [Note that while the actual point group at U sites in 
UPts is -D3/1, the allowed crystal field spectrum is essentially identical to that of Dq^.] 

A single-site Hamiltonian H is defined coupling the quadrupolar operators with an 
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electric-field gradient V-E and coupling the magnetic moment with an external magnetic-field 
H, at the uranium site: 

H = Ho-H-M- Y.i^c.Ep)Qf,^ , (5) 

where Hq is the Hamiltonian for the uranium ion in the presence of the crystal field only 
(see Fig. |1|): 

Ho = A\B){B\. (6) 

We are measuring energies from the doublet ground-state \E±). 
For in-plane H and E, the Hamiltonian Eq. (^) is written as 

H = Hq + AJ+ + BJl + h.c. , (7) 

where 

A = -^iH,-tHy) (8) 

and 

B = -^ (dxxE - dyyE - 2ld^yE) . (Q) 

The diagonalization of this Hamiltonian shows that there is a crossing of levels depending 
upon H and E and an interplay between the field gradient and magnetic field which depends 
crucially on their orientation with respect to one another. In the model defined in the next 
section, the coefficients A and B are interpreted as molecular fields acting on the uranium 
sites. Since those molecular fields depend on the magnetic and quadrupolar moments, a 
self-consistent calculation is required to obtain these moments at each temperature. 

B. The Intersite-coupling Hamiltonian 

In order to determine the molecular fields H and E present in Eq. (|^) we have used 
a model Hamiltonian which is Heisenberg like in the magnetic moments as well in the 
quadrupolar moments, defined by 



<i,j> 

+ ifTrQ, ■ Q, + IlTrQ, ■ Q,) , (10) 

where I'^ (/^) is the couphng between the nearest-neighbour magnetic moments in the 
basal plane (c-axis direction) of the hexagonal lattice; (J^) has the same meaning for the 
quadrupolar moments. Tr means trace over the product of the quadrupolar tensors. 

We are interested in studying the antiferromagnetic order of the model defined by 
Eq. (p!0D, so that the couplings I'^ and are negative. Only ferroquadrupolar coupling Jq 
has been considered. We have examined the effect of a coupling of the form TrMi ■ Qj ■ 
in the Hamiltonian Eq. (plQl) . It does not change qualitatively the mean-field phase diagram 
of the model, since the couplings already present in Eq. ([I0|) give rise to a term proportional 
to M'^Q in the free energy (see Eq. (|1^)). Hence we will drop this coupling from the Hamil- 
tonian. The calculations are restricted to the case when M lies in x-direction (that is the 
case for UPts, as showed by neutron scattering data [§]), although it is possible to consider 
the general case. 



III. MEAN-FIELD THEORY 



A. Mean-Field Hamiltonian 



The mean field approach consists in replacing the operators in the Hamiltonian by their 
mean values plus the respective fluctuations, and keeping the fluctuations only up to first 
order. Denoting the mean value of the magnetic and quadrupolar moment by M = gjfiBd^ 
and Q = bq, respectively, where the matrix elements a and b are defined in Eq. (^ and 
Eq. (H), the mean field apphed to Eq. (p!OD results in 

I = ^ - N{Am + Bq) + 2 ^(Amf + BgH , (11) 

i 

where N is the number of sites of the hexagonal lattice, 

A = ilM±jMKm ^ J^m , (12) 
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and 



B = ^ -6J,<i . (13) 



The last equalities, Eq. ([1^) and (|13]), define the dimensionless magnetic coupling Jm and 
the dimensionless quadrupolar coupling Jg, respectively. We are using deliberately the same 
notation A and B here as in Eq. (Rf) and (01) since in both cases the meaning of those 



parameters are the same, although in Eq. (|T2D and (|l^) A and B depend on the order 
parameters rh and q. 

The site-independent terms of H are called the condensed energy; the other terms define 
a sum over an effective one-site Hamiltonian Hi, which can be written in terms of the 
components of the total angular momentum as 

= Ho^ + A{Ji+ + J,_) + B{Jl + Jl) . (14) 

Diagonalizing Hi in the subspace defined by the electronic levels shown in Fig. (|I]) we 
find the following expression for the free-energy per site and in units of A: 

/ = -Jmm^ + QJqq^ - nYiie^P^^"^" + 2e-^^^-^''^'^'^ cosh/5C) + t ln(2 + e"'') , (15) 

where (3 = 1/tis the inverse of the dimensionless temperature t = T/A. C is given by 



C=^{^-±^y + 2Jlm^. (16) 

From the free-energy, we obtain the following coupled equations for the order parameters 
rh and q : 

^ ^ _^JmJn sinh/jC 

^~ C exp[l3{l - 18Jgg)/2] + 2cosh/5C ' ^ ^ 

and 

lcosh/5C+ (l + 6Jgg)/(2C)sinh/5C-ea;p[/?(l - 18J5g)/2] 
^~2 exp[(5{l - 18Jgg)/2] + 2cosh/5C ' ^ ^ 

Before solving equations ( p!7| ) and (|TB|) numerically, we describe a few limits where it is 
possible to study these equations analytically. There is no solution other than the trivial 
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m = and g = if Jm < — 1/2 and Jq = 0. A pure magnetic phase (m 7^ 0, g = 0) 
does not exist, but a pure quadrupolar phase (m = 0, g 7^ 0) exists for | J^l < (1 + 3Jg)/2 
and Jq 7^ 0. In that phase, for T << A one has q ~ 0.5 tanh (6/3Jgg), which imphes the 
quadrupolar degrees of freedom are playing the role of a pseudospin 1/2. In this regime, a 
second-order critical line appears at tcq ~ SJq. The other regimes of Equations (|T^ and (|TBp 
are studied numerically. The next section discusses the phase diagram obtained by studying 
the stabilities of the solutions rh and q for different temperatures and values of the coupling 
coefficients Jm and Jq. 

B. Phase Diagrams 

Fig. shows the phase diagram of the model defined by Eq. ([TT|). The main features of 
this diagram are the enhancement of the critical temperature and the interchange between 
first and second-order transitions when the quadrupolar coupling Jq is introduced. 

We first consider the case when Jq = 0. Both order parameters, rh and q, are critical 
simultaneously and despite the fact that we do not have quadrupolar coupling in this case, 
the ground state described by the Hamiltonian H, Eq. (0), has quadrupolar degrees of 
freedom and for different temperatures the stability of this ground state against the intersite 
magnetic coupling causes the tricritical point (denoted by TP), separating the first-order 
and second-order transitions. The coordinates of this tricritical point are Jm = —1.05 and 
tc = 0.300. The line of second order is obtained from the limit rh and q going to zero in 
Eq. ( p!?!) and Eq. (p^) and is given by: 6^^= = (Jm — 1/2)/ {Jm + 1)- Near this line rh follows 
the usual (tc — t)^/^ behavior characteristic of the order parameter in mean-field theory; q is 
proportional to m^, which is characteristic of the secondary order-parameter. The first-order 
line extends from the tricritical point to Jm = —1/2; below this point no order is found for 
either rh or q. 

For non-zero Jq we have simultaneously first-order transition for both order parameters 
only above the critical-end-point, denoted CEP in Fig. ^ Below this point the transitions are 
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separated in one line of second-order, purely quadrupolar transition at t^q = 6Jq/ (2+e~^/*'=9), 
and another line of first-order purely magnetic transition. This latter line has a tricritical 
point, denoted TP, becoming second-order below this point and finishing at t = and 
Jm = — (1 + ^Jq)/'^] below this value of Jm one has only quadrupolar order. In Table | 
we have listed positions of the points CEP and TP for a few values of Jq. Fig. ^ shows 
typical behavior of both order-parameters with temperature, in different parts of the phase 
diagram. 

Besides changing the values of the critical temperatures and separating the phases below 
the CEP, the intersite quadrupolar coupling has a more dramatic effect on the phase dia- 
gram, inverting the nature of the phases with respect to the TP: the first-order transition 
occurs below the TP when Jq = 0, and exactly the opposite occurs when Jg ^ 0. We be- 
lieve this complicated effect of the indirect magnetic-quadrupolar coupling generated by the 
Hamiltonian Eq. (0) exists even for a very small quadrupolar coupling. This is confirmed 
by Landau's mean field theory, which is discussed in the Section 



C. Specific-Heat Jumps 



Table |IT| shows numerical calculations of the value of the specific heat jumps, AC, for 
the second-order magnetic transitions for different quadrupolar coupling Jq. For zero J^, AC 
is never less than QJ/mol — K. For non-zero Jg, AC is a very sensitive function of Jm-, giving 
values of AC around 2J/mol — K for Jq = 0.05, and Jm = —0.6, which is the experimental 
value of AC for the compound U1_1.Th1.Pt3. For these Jm and Jq, and considering that 
the Neel temperature is equal to 6K for this compound, the model also gives the value of 
energy splitting A of the crystal field as 70K, which is the order of magnitude expected for 
this material. Naturally we will have corrections to the specific-heat jump when including 
explicitly the interaction with the conduction-band electrons. If these corrections are not 
too large, the model shows that a very small quadrupolar coupling, about one order of 
magnitude smaller than the magnetic coupling, is necessary to obtain the magnitude of the 
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specific-heat jump. 



IV. LANDAU'S MEAN-FIELD THEORY 

Contrary to the previous section, which treated the molecular-field-theory free energy 
exactly, now we want to use the Landau's mean field theory, which is valid in the limit of 
small order parameters, i.e., near a second-order transition. This theory, although limited, 
can describe some of the features of the phase diagrams. 

A. Free Energy 

The Landau's mean-field theory uses an approximation for the free energy, valid near 
the region where the order parameters m and q are small. Expanding the free energy given 
by Eq. (|15D around m = and g = we have: 

/ = aynfa^ + f^mrfi^ + Imrn^ + oiqt + f^q^ + Smqrh'^q + ernqfa^t + ••• , (19) 

where the coefficients are functions of the temperature t and of the couplings Jm and Jq. 
They are listed in Table |T[ 

For Jq = 0, only the three first terms on the right hand side of Eq. (0) survive 
and the critical temperatures for second-order transitions are obtained from the equation 
o^mitc, Jm) = 0, whose solutions agree with the numerical results. The coefficient Pmit, Jm) 
is negative (positive) for temperature less (great) than 0.30032, and 7^ is positive (negative) 
for temperature less (great) than 0.46370, for all Jm- That means we have first (second) order 
for critical temperatures less (great) than 0.30032. The point where both coefficients, am and 
Prn, vanish simultaneously determines the tricritical point (TP) and that is at txp = 0.30032 
and Jm^ = —1.07140, also in agreement with the phase diagram Fig. (^. 

The quadrupolar-coupling Jg complicates the analysis of the free energy Eq. (|1^). Now, 
the second-order magnetic transition (dashed line in Fig. 0) has non-zero quadrupolar order- 
parameter q = qo, so the free energy in fact should be written as 
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/ = a'^{t, Jm, qo)^^ + P'mit, Jm, qo)m'^ + ••• , (20) 

where go is the equihbrium value of q on the critical line (m = 0), and is given by go = 
(— ag/2/5g)^/^. The critical temperature for the magnetic transition now is the solution of 
a'^itc, Jm,qo) ~ ttm + iqo + ^qi + ••• = 0, which gives tc close to the numerical results, 
mainly near (but at left of) the tricritical point where go is not to big. To calculate (3'^ 
needs some care. Using the expansion (^) (3'^ will have divergent terms at t^p, which 
invalidate its expansion around the tricritical point since this point occurs always near t^p. 
For Jq = those terms do not exist and = Pm- We have computed P'^ through its 
definition: = {d'^ f / dm^) / 41 , calculated at m = and g = go, where go is obtained from 
the numerical solution of the Eq. (\[7\} and Eq. ([18|) . showed no divergence and, although 
Prn can be negative, is positive for t < t^p, which explains why we have a second-order 
transition when Jg intead of a first-order one when Jg = 0. Unfortunately, Landau's 
mean-field theory does not work well right above the tricritical point TP since m does not go 
to zero, and in fact rh has a big jump at the critical temperature and so g has a discontinuity 
at this line (see Fig. ^). For large Jm, however, the jumps in both order parameters at the 
critical line become small, thus the expansion (p!9| ) gives a reasonable description of the 
model in this limit. 

V. SUMMARY AND CONCLUSIONS 

We have calculated the phase diagram of a simple model which includes intersite 
quadrupolar coupling between uranium ions in the hexagonal lattice. The model has some 
of the physics found in the hexagonal uranium-based materials, like the in-plane van Vleck 
magnetism, magnetic character along the c-axis and degrees of freedom to allow the possi- 
bility of quadrupolar-Kondo effect. We have shown that the intersite quadrupolar coupling 
between uranium ions enhances the critical temperature and plays the fundamental role in 
defining the nature of the transitions in the phase diagram. 

Based on the magnitudes of the specific-heat jumps and the magnetic moment fh given 
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by the model around the tricritical point (TP) when Jq is of order |Jm|/10, we speculate 
this region as most probable to find the hexagonal uranium-doped compounds. We are 
performing a model calculation based on the Coqblin-Schrieffer [|I0| formalism to determine 
whether this coupling constant ratio can be realized in practice. 

Naturally, the model defined here represents just the first step in the development of a 
model for the antiferromagnetic transition in U-heavy-fermion materials. We should next 
introduce explicitly the conduction electrons, in hence the Kondo effect. The natural ques- 
tion is what features of the mean-field phase diagram will survive. The structure of the 
phase diagram near the tricritical point for non-zero Jq appears essentially due to the terms 
5mqfn?<i and emq'^n^'f in the free-energy, Eq. ([T5|), generated by the Hamiltonian (pTj). Gu- 
fan et al. using Landau's mean-field theory have analysed the phase diagram of a free 
energy of the form given by (0), and they also found the same topology we have around 
the multicritical points. Those terms in our model are always present as long as we have 
ferroquadrupolar coupling and the hexagonal symmetry. 

The Kondo effect will not alter the structure of the Ginsburg-Landau free energy which is 
dictated purely by symmetry considerations. However, will renormalize the GL coefficients. 
We don't expect the phase diagram to be strongly modified in the absence of intersite 
quadrupolar coupling, since the physics of the magnetic susceptibility remains that of van 
Vleck with only slight modification of the low temperature values relative to the case of 
zero Kondo temperature. However, once we turn on intersite quadrupolar coupling, we will 
dramatically renormalize the coefficients. In particular, wherever a 1/T appears, in our 
expressions, we reahze that at a first pass this is the single-site quadrupolar susceptibihty 
of the U ions, and given the quadrupolar Kondo effect, which is two-channel, this will 
diverge logarithmically |T^. Using the single site susceptibility of the two-channel model as 
a guide, we will effect the replacements 5mq ~ 1/T ^ 1/Tk\^{Tk /T) and e^q ~ 1/T — >• 
1 /Tft- In (Tft:/T) where Tk is the Kondo temperature. This will not remove the second- 
order quadrupolar and magnetic transitions in our phase diagram, but will renormalize 
the transition temperatures downwards. As a result of this we anticipate that we may 
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obtain similar reduced moment values and specific heat jumps to UPts with larger intersite 
quadrupolar couplings. 

In addition to this exploration of how the Kondo effect modifies our already interesting 
results, we would like to explore the effects of magnetic field through magnetostrictive cou- 
pling to see whether we can understand the metamagnetic transition observed in pure UPta 
below 20K in fields of order 20T [O]. 
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TABLES 

TABLE L Coordinates of the tricritical points (TP) and critical end points (CEP) in the phase 
diagram of the Fig. ^, for different quadrupolar coupling Jq. 



Jq tTP -Jm^ tcEP -J^^^ 

0. 0.30 1.06 — — 

0.05 0.12 0.64 0.14 0.70 

0.1 0.27 0.79 0.29 0.83 

0.2 0.52 1.01 0.55 1.10 



TABLE II. Coefficients of the expansion of the free energy defined in Eq. ([T9|). Here 

am = —^dJm sinh/?/2 — 

= 2d^Ji[{2 - (3)eP - 5/? - 1 - e-^] 

= 8J^{3(3^d^ - M) sinh[/3/2] + 8J^pd^{5 + e^) 
ag = QdJqef^/^{2 + e"^ - 6Jqf3) 
(3q = 108^2 Jg4/?3(4e^ - 1) 
5mg = l2dJlJq{(5e^l^ - 2sinh/3/2) 
eraq = 36^^ (4/?e'3/2 - 46^3/2 _ 3/32 ^ 2/3 + 2 + 2e'^) 
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TABLE III. Specific-heat jumps AC of the second-order magnetic transitions on the phase 
diagram Fig. ^, for a few values of the quadrupolar coupling Jq. The magnetic coupling are 
near the minimum value to give second-order transitions. 



0.00 1.10 6.52 

0.05 0.60 1.50 

0.1 0.68 0.92 

0.2 0.85 1.31 
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FIGURES 

FIG. 1. Level structure of the model considered in this work. The doublet ground-state 
\E±) = ti| lb 4) + ti| =F 2) is magnetic coupled to a singlet excited-state V^\B) = |3) -|- | — 3), which 
has energy A. Here \m) is the eigenstate of the z-component of the total angular momentum J = 4. 
The other D^^ crystal-field levels are not included, thus restricting the model to temperatures less 
than A. 



FIG. 2. Mean-field phase diagram of the model defined by the Hamiltonian Eq. (11). Jm and 
Jq are the dimensionless magnetic and quadrupolar coupling between the nearest-neighbouring 
sites, respectively, and tc = ksTc/A is the dimensionless critical temperature. At right of the 
critical end point (denoted by CEP), both order parameters have the same critical temperatures. 
At left, the transitions are separated into one pure second-order quadrupolar (horizontal straight 
line) and one pure first-order magnetic transition (dashed line). This last one ends at the tricritical 
point (denoted TP) where a line of second-order magnetic transition begins and ends at i = and 
Jm = —0.5(1 -|- 3Jq). The quadrupolar coupling enhances the critical temperature tc, increases 
the absolute value of the minimum magnetic coupling Jm necessary to have magnetic order, and 
reverses the nature of the phase transitions respect to the tricritical point. 

FIG. 3. Temperature dependence of the magnetic moment m and quadrupolar moment q, 
obtained from Eq. (17) and Eq. (|l^). The temperatures on the x-axis are measured in units of 



the crystal-field splitting A, and are not in the same scale; the units of all other quantities are 
defined in the text. For Jq = 0, Fig. (a), (b) and (c) q has the same critical temperatures as fh, 
and behaves like a secondary order-parameter. For Jq = 0.1, Fig. (d), (e) and (f) the quadrupolar 
transitions are separated from the magnetic ones, except for \Jm\ > \Jm^^ = 0.83. For the same 
value of Jm the nature of the phase transitions changes according to whether Jq is present or not. 
No magnetic order is found without quadrupolar order. 
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